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Abstract 
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•7- . We derive a formula that simplifies the original asymptotic iteration method formulation to find 

00 ' the energy eigenvalues for the analytically solvable cases. We then show that there is a connection 

between the asymptotic iteration and the Nikiforov-Uvarov methods, which both solve the second 

O ■ order linear ordinary differential equations analytically. 
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Analytical solutions of the radial Schrodinger equation are of high importance in quan- 
tum mechanics since the wave function contains all the necessary information to describe a 
quantum system fully^^"^^'. The asymptotic iteration method (AIM)[^'^1 and the Nikiforov- 
Uvarov method^^^ are two different approaches to solve the resulting Schrodinger equation 
and they have been successfully applied to a wide variety of problems f^^^^'. 

In this work, we derive a formula which simplifies the original AIM formulation to find 
the energy eigenvalues for the analytically solvable cases. This formulation shows that there 
is a connection between the AIM and Nikiforov-Uvarov methods regarding the solution of 
differential equations. 

AIM'^1 is proposed to solve the second-order differential equations in the form 



v" = Xoir)y' + So{r)y, 



(1) 



where Ao(r) 7^ and So(r), Ao(r) are in Coo(a,b). Eq. (1) as given in Ref. [1] has a general 
solution 



y(r) = exp 



a{r )dr 



//,,,// \ , / 



C2 + C1 / exp / (Ao(r ) + 2a{r ))dr dr 



(2) 



if A; > 0, for sufficiently large k, a{r) values can be obtained by 



Afc(r) Afc-i(r) 



a r 



(3) 



where 



Afc(r) = A^_i(r) + Sfe_i(r) + Ao(r)Afc-i(r), 

Skij-) = s'f._^{r) + So{r)\k-i{r), k = 1,2,3, 



(4) 



In this method, the idea is to convert the radial Schrodinger equation into the form of 
Eq. (1) for a given potential. Then, So(r) and Ao(r) are determined, and Sfc(r) and Afc(r) 
parameters are calculated. The energy eigenvalues are determined by the quantization con- 
dition given by the equation: 



Sk{r) = Afc(r)sfc_i(r) - Afc_i(r)sfe(r) = 0, A; = 1, 2, 3, ■ ■ • . 
If we rearrange Eq. (1) by introducing new forms for Ao(r) and So(r) as follows: 



Ao(r) 



Tir)_ 



so{r) 



In 

air) 



(5) 



(6) 



where r and a are new functions and 7„ is a constant with which we quantize the energy. 
Inserting them into Eq. (1), the second-order differential equation takes the form 

V = — rrV rrV- iV 

a[r) (j(r) 

An essential ingredient of AIM is the quantization condition given by Eq. (5). Every 
particular value in Eq. (5) corresponds to an energy eigenvalue. Therefore, if we first take 
k = 1 and n = 0, put them into Eq. (6) and then use Eq. (4) for 70 7^ 0, the quantization 
condition given by Eq. (5) reduces to 

r'(r) + 7o = 0, (8) 

which should hold for all r. This means that T{r) is a polynomial of the first degree at most. 
If we choose now n = 1 and k = 2, using Eq. (4), the quantization condition given by Eq. 
(5) reduces to 

7i + r'(r) {-fi + 2T'{r)+a"{r))-T"{r){T{r) + a'{r)) =0, (9) 

which should also hold for all r. Since we have already established that T{r) is a polynomial 
of the first degree at most, thus, T"{r) = 0. Assuming the spectrum of our problem to be 
non-degenerate, i.e. 71 7^ 70, and taking into account Eq. (8), we conclude that a'"{r) is a 
polynomial of the second degree at most. In such a case, Eq. (7) and, hence, Eq. (1) are of 
hypergeometric type (We point out that the comments in this part have been made by an 
unanimous referee of the Phys. Lett. A for another paper of ours). 

Already knowing the properties of T{r) and cr(r) and by using AIM procedure given 
by Eq. (4), we can calculate A„(r) and s„(r) as follows, in order to formulate the energy 
eigenvalue equation, 

7(r(r) + a^(r)) „ ,,^ _ r(r) {rir) + g^r)) - ajr) (7 + T'ir)) 
M{r - -—^ , Si[r — -—T , (^iUj 

■ ■ ■ etc 
If we use the termination condition of AIM given by Eq. (5), the energy eigenvalues are 



obtained as follows 



So 
Ao 


Sl 

Ai 


Sl 

Ai 


S2 
A2 


S2 
A2 


S3 
A3 


S3 
A3 


S4 
A4 


etc 





70 = (11) 

^ = f ^ 71 = -r'(r) (12) 

72 = -2r'(r) - a"(r) (13) 

73 = -3T'(r)-3a"(r), (14) 



which can be generalized as 



7„ = -nr'ir) - !^i!^a"(r) (15) 



This formula gives directly the energy eigenvalues for differential equation (7). It replaces the 
iterative calculations needed to find the energy eigenvalues in the original AIM formulation. 
Also, it is exactly the same as the Nikiforov-Uvarov energy eigenvalue equation. This is 
briefly derived as follows: 

The Nikiforov-Uvarov equation is written as 

V^(r)" + ^^\r) + ^Mr) = 0, (16) 

where a{r) and d'{r) are the polynomials of the second-degree at most, and f(r) is a 
first-degree polynomial '^'"^'^^^l. Hence, from Eq. (16), the Schrodinger equation or the 
Schrodinger-like equations can be solved analytically by this method. In order to find a 
particular solution of Eq. (16), the following transformation is used: 

^(r) = (p{r)y{r). (17) 

It reduces Eq. (16) to an equation of hypergeometric type, 

a{r)y" + T{r)y' + 7?/ = 0, (18) 

and (f){r) is defined as a logarithmic derivative in the following form and its solutions can be 
obtained from 

(P'{r)/(P{r) = n{r)/a{r) (19) 

The function n and the parameter A required for this method are defined as follows: 




^W = ^^±\/^^ -a + A;a, (20) 



A = A; + tt'. 



(21) 



On the other hand, in order to find the value of k, the expression under the square root must 
be the square of a polynomiaL Thus, the eigenvalue equation for the Schrodinger equation 

becomes 

, n{n - 1) 
An = —nr a . [22) 



Next, we present some applications. For classical differential equations, in order to show 
the applicability of the above formula, let us consider the classical differential equations 
which give polynomial solutionst^^' . The results are presented in Table 1. 

For some other applications reported previously in the literature by using AIM, we may 
take the following examples: 

In the Morse potential case, if we take the transformed equation (31) of Ref. [16] as 

dVn(r) f2i3y-2e-a\dfnir) , f2eP + aP-2(3^- 



dr"^ 



ar 



dr 



+ 



ra^ 



Ur) 



where 



— rij-) = 2(3r — 2e — a, a = ar, — 7„ 
Using Eq. (15), the eigenvalues turn out to be 



2ep + ap- 2(3"' 



a 



^n = l3-{n+-)a n = 0,1,2,3, 



(23) 



(24) 



(25) 



For the deformed Hiilten potential case, if we take the transformed equation, Eq. (18) of 
Ref. [16] can be rewritten as 



d'f{r) 
dr'^ 



{2e + 3)qr - {2e + 1] 
r(l — qr) 



df{r) 
dr 



2eq + q- P^ 
r(l — qr) 



fir). 



If we compare this equation with Eq. (6), we obtain the functions 



r(r) = (2£ + 3)gr- (25+ 1), a{r) = -r{l - qr), -jn = '2eq + q-p'^ 



(26) 



(27) 



Using Eq. (15), the eigenvalues turn out to be 

1 ^g(n+ 1)2-/52 



n = 0,l,2,3. 



(28) 



(2n + 2) V q 

In the Kratzer potential case, if we take the transformed equation, Eq. (14) of Ref. [17] 
can be rewritten as 



/•(,) . 2 (. - (^) /» + (?(^±i^i^) /(,■). (29) 

If we compare this equation with Eq. (6), we obtain the functions 

- T{r) = 2 (er - (A + 1)) , a = r, --fn = 2(A + 1)£ - A. (30) 

Using Eq. (15), the eigenvalues turn out to be 

'g(n + 1)2-/32 



(2n + 2)V 9 ;• " = °'l-2-3.- (31) 

Next, we discuss the eigenfunctions. If we insert Eq. (15) into Eq. (7) and use Eqs. (2) 
and (3), we can obtain the wave functions as follows: 

Voir) = 1, (32) 

yi{r) = r, (33) 

y^[r) = T^ + Ta' + T'a + aa", (34) 

y^(r) = r^ + 3rV + 2rcT'2 + 3rr'a + 4rW + 5raa" + 6aaV, (35) 



These iterative functions can be generalized in a closed form as follows: 

1 rf" 
Vnir) = ^^ HrTpir)] , n = 0, 1, 2, 3, ■ ■ ■ , (36) 

where p{r) has to satisfyt^' 

[a{r)p{r)]' = r{r)p{r). (37) 

If we apply Eq. (36) to the deformed Hulten potential case as an example, we obtain 



Voir) = 1, (38) 

l/i(r) = (2e, + 3)qr-{2e, + l), (39) 

2/2(r) = 2 [{2€l + 352) + (£2 + 2)(2£2 + 5)gV2 _ ^^^ + i^^^^ + 2)qr + l] , (40) 

ysir) = 2(£3 + 3) (253 + 5) (2^3 + 7)gV - 6(^3 + 3) (2^3 + 3) (2^3 + 5)qh^ + (41) 
+2(3£3(4£3(^3 + 5) + 31) + 45)gr - 2(£3(4£3(^3 + 3) + 11) + 3), 



These equations can be written in a generalized way as a hypergeometric function, 

yn{r) = (-l)" ^^2et+y^ ^'^^ ^~''' ^^^" + ^ + 2) 5 {2^" + U; ^0 . (42) 

If Eq. (42) is compared with the resuhs of Ref. [16], it can be seen that they are exactly 
the same eigenf unctions. 

In summary, we have derived a formula for the analytically solvable potentials, which is 
given in Eq. (15). This formula replaces the iterative calculations and reduces the calculation 
workload in the original AIM formulation to find the energy eigenvalues. We have also shown 
that there is a close connection between the AIM and the Nikiforov-Uvarov methods to solve 
differential equations. 
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Differential equation 


— T 


a 


-7n 


Eigenvalue 


Cauchy-Euler^ 


a{r — b) 


(r - af 


/? 


/3 = n{n — 1 — a) 


Hermite^" 


2r 


1 


-2A; 


k = n 


Hermite^'' 


ar + b 


1 


c 


c = —na 


Laguerre 


r-l 


r 


a 


a = —n 


Confluent^ 


br — c 


r 


a 


a = —nb 


Hypergeometric 


{a + b+l)r -c 


r(l-r) 


ab 


a = —n or b = —n 


Legendre 


-2r 


r2-l 


m{m + 1) 


m = n 


Jacobi 


{a + p + 2)r + P + a 


l-r2 


-7 


7 = n{n + a + P + 1) 


Cliebyshev^'* 


r 


l-r2 


— m 


m = V? 


Cliebyshev^^ 


3r 


l-r2 


— m 


m = n{n + 2) 


Gegenbauer 


(l + 2fe)r 


(l-r^) 


-A 


A = n(n + 2A;) 


Hyperspherical 


2{1 + k)r 


(l-r2) 


-A 


X = n{n + l + 2k) 


Bessel^" 


-2(r + l) 


r2 


7 


7 = n{n + 1) 


Generalized Bessel^ 


-{ar + b) 


^2 


7 


7 = n{n + a — 1) 



TABLE I: Application of Eq. (15) to classical differential equations which give polynomial solutions. 



